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ABSTRACT 


Stability  map  for  a  gas-bearing  supported  rotor  can  be  constructed  from  the 
periodic  dynamic  perturbation  solution  of  the  gas  bearing  equation  about  the 
equilibrium  condition.  A  complete  frequency  range  of  the  perturbation  solution 
is  required.  Effective  system  stiffness  and  damping  can  be  defined  in  terms  of 
the  perturbation  solution  of  the  gas  bearing.  The  condition  of  neutral  stability 
corresponds  to  the  vanishing  of  the  effective  damping,  yielding  the  characteristic 
frequency,  and  resonance  at  the  characteristic  frequency  with  the  equivalent 
stiffness,  yielding  the  critical  mass.  An  excess  of  rotor  mass  causes  instability 
if  the  effective  damping  Increases  with  frequency,  at  the  characteristic  frequency, 
and  conversely.  For  unloaded.  Journal  bearings  possessing  rotation  L  s'-mmetry, 
the  stability  map  can  be  constructed  from  their  steady-whirl  solutions.  .\n  example 
illustrating  the  latter  case  is  given  in  terms  of  a  herringbone -grooved  Journal 
bearing.  Ihis  method  of  stability  analysis  is  applicable  to  both  thrust  and 
Journal  bearings  for  both  whirl  and  pneumatic-haomer  instabilities. 


INTRODUCTION 


Certain  types  of  fluid  film  bearings  are  capable  of  causing  some  form  of  dynamic 
Instability  of  the  bearing-rotor  system.  The  best  known  example  ds  the  whirl 
Instability  of  self-acting  Journal  bearing,  which  first  became  understood  for 
oil-lubricated  journal  bearings  in  Poritsky's  "Theory  of  Oil  Whip"  (Ref.  1). 

The  same  phenomenon  also  exists  with  self-acting  gas-lubricated  Journal  bearings 
and  has  been  referred  to  as  "ha If- frequency  whirl"  or  "fractional- frequency 
whirl".  Although  experimental  data  on  the  whirl  instability  of  rotors  supported 
in  gas  Journal  bearing  (Refs.  2,  3,  4,  and  5)  are  quite  similar  to  those  of 
rotors  supported  in  oil  Journal  bearings,  attempts  to  extend  the  "oil-whip  theory 
to  gas-bearing  systems  encountered  considerable  difficulty.  This  is  primarily 
because  the  gas- lubrication  equation  explicitly  contains  the  time  derivative  of 
gas  density  (Ref.  6).  While  an  analogy  of  this  equation  with  the  heat  diffusion 
problem  is  quite  apparent  (Refs.  7  and  8),  the  significance  of  this  term  in  the 
dynamics  of  rotor-bearing  systems  has  been  an  object  of  controversy  for  some  time 

From  the  standpoint  of  obtaining  the  complete  mathematical  solution  of  the  gas- 
lubrlcatlon  equation,  it  is  beyond  doubt  that  an  initial  condition  should  be 
specified.  On  the  otherhand,  it  is  intuitively  plausible  that  the  stability  of 
the  rotor-bearing  system  should  not  have  a  history  dependence.  Earlier  works  on 
the  subject  reflected  the  latter  view  and  neglected  the  time  dependent  term  from 
the  gas-lubricated  equation  altbgether  (Refs.  9  and  10);  clearly,  the  accuracy 
of  these  results  are  not  a  priori  known.  Subsequently,  it  has  been  shown  that 
time  dependence  of  the  gas-film  pressure  can  be  separated  into  an  "initial 
transient",  which  is  directly  related  to  the  specific  initial  condition  and  would 
always  attenuate,  and  an  "implicit  time-dependent"  tera,  which  interacts  with  the 
motion  of  the  rotor  and  controls  the  dynamic  stability  of  the  system  (Ref.  8). 
Also,  there  is  a  formal  method  to  uncouple  the  time -dependent  gas-luhrication 
equation  from  the  dynamical  equation  of  the  rotor  (Ref.  8).  This  method  of 
uncoupling  is  equivalent  to  a  Taylor  series  expansion  in  the  time  coordinate 
and,  therefore,  is  subject  to  truncation  error  which  is  not  readily  controllable 
and  can  significantly  affect  the  stability  analysis  (Ref.  11).  Hore  recently 
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the  stability  of  a  rigid  rotor  in  plain,  cylindrical,  self-acting,  gas-lubricated, 
journal  bearings  have  been  studied  quite  thoroughly  with  the  time-dependent  effects 
properly  considered  (Refs.  12,  13  and  lA).  Here,  peculiarities  in  the  stability 
analysis  of  gas-bearing  systems  are  again  manifested.  Where  the  method  of  Galerkin 
is  utilized,  the  characteristic  equation  of  the  system  (Refs.  12  and  14)  is  a 
polynominal,  and  Routh-Hnrnitz  criteria  can  be  used  to  determine  the  condition 
of  stability.  Since  the  accuracy  of  the  method  of  Galerkin  is  improved  by  an 
increased  number  of  terms  used,  the  corresponding  characteristic  equation  also 
becomes  of  a  higher  order.  Where  the  finite-difference  method  is  employed  (Ref. 
13),  the  state  of  periodic  motion,  which  borders  the  domain  of  unstable  operation, 
can  be  found.  In  this  analysis,  one  does  not  know,  except  with  the  additional 
effort  of  orbit  calculation  whether  one  or  the  other  side  the  neutral  state  is 
stable. 

Now,  stability  problems  are  not  limited  to  the  self-acting  plain  Journal  bearings. 
The  externally-pressurized  gas- lubricated  thrust  bearing  can  exhibit  a  form  of 
self-excited  oscillation  commonly  called  "pneumatic  hammer"  (Ref.  15),  which  Is 
also  possible  for  an  externally-pressurized  Journal  bearing.  External- 
pressurization  can  reduce  but  will  not  eliminate  the  hazard  of  whirl- 
instability  of  high-speed  Journal  bearings  (Refs.  16,  17  and  18).  A  variation 
of  "pneumatic  hammer"  has  been  discovered  for  the  spiral-grooved,  self-acting, 
gas- lubricated  thrust  bearing  (Ref.  19).  Many  'Vthirl-free"  Journal  bearing 
configurations  have  been  proposed  (Refs.  20,  21,  22,  and  23). 

The  importance  of  stability  analysis  in  the  application  of  gas  bearings  has 
already  attracted  attention.  However,  a  general  method  for  performing  the 
stability  analysis,  %ihich  is  convenient  to  use,  suggestive  in  physical  meaning, 
free  from  arbitrary  assumptions,  and  imdependent  of  the  particulars  ot  the 
bearing  configuration,  is  still  missing  from  public  literature.  The  present 
work  is  intended  to  fill  this  gap,  at  least  in  part. 

It  will  be  shown,  that  the  periodic  dynamic -perturbation  forces  of  the  gas  bearing 
can  be  used  to  construct  the  effective  damping  and  stiffi^ss  coefficients  of  the 
dynamical  system,  that  the  state  of  neutral  stability  occurs  at  a  characteristic 
frequency  which  is  determined  by  requiring  the  effective  damping  coefficient  to 
vanish  and  when  the  rotor  smss  has  a  critical  value  which  is  in  resonance  with 
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the  effective  stiffness  at  the  characteristic  frequency,  and  that  an  increment  of 
the  rotor  mass  beyond  the  critical  value  will  cause  instability  if  and  only  if 
the  effective  damping  coefficient  increases  with  frequency  (at  the  characteristic 
frequency).  Thus,  the  key  to  this  method  of  stability  analysis  is  whe  knowledge 
of  the  bearing  reactions  to  perturbed  periodic  motions  for  the  entire  frequency 
spectrum. 

Derivations  will  be  carried  out  for  a  single  mass  supported  in  a  gas  bearing 
having  two  degrees  of  freedom.  Simplifications  for  the  special  cases  of  a  single 
degree  of  freedom  bearing  (e.g.,  thrust  bearing)  and  an  isotropic  bearing  (e.g., 
unloaded  Journal  bearing)  will  also  be  given.  To  illustrate  the  application  of 
this  method,  the  stability  map  of  an  unloaded  herring-bone  grooved-journal  bearing 
will  be  constructed. 
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A  GENERAL  THEOREM  FOR  THE  STABILITY  OF  A  RIGID  ROTOR  IN  FIXED  JCHIRNAL  BEARING 


The  gas  lubrication  equation  for  a  given  Journal  bearing  has  a  solution  dependent  on 

the  motion  of  the  journal.  Consider  Fig.  1.  Let  the  Instantaneous  displacement  of 

the  Journal  from  Its  equilibrium  position  be  (e  »  Assume 

X  y 

e  ^  e  ^  <  <  (1) 

X  y 

C  <'x-  ‘y>  ■  <2) 

where , 

T  ■  cot 

s  »  \  +  iv  (3) 


Then  the  lubrication  equation  can  be  llnearlced  and  yield,  as  final  results,  the 
perturbed  dynamic  reaction  of  the  bearing  to  the  Journal  motion.  It  is  customary 
to  make  the  bearing  reaction  non>dimenaional  with  the  normalixing  factor  p^LD. 
Linearisation  leads  to  the  useful  simplification  that  the  perturbed  bearing 
reaction  is  directly  proportional  to  the  motion.  Thus,  we  can  write  In  Argand 
notation  the  matrix  equation  * 

f  -  2(s)€  (^) 


where,  f  Is  the  column  vector 


^e  homogeneous  fora  of  Eq.  (4)  impliea  that  the  Initial  transient  of  the  fluid 
film  pressure  has  been  n«.glected.  For  justification  of  this  aasumption,  see 
Ref.  8. 


5- 


which  is  related  to  the  actual  components  of  the  perturbed  bearing  reaction  according 
to 


_1 _ 

p  LD 


Z(s)  is  a  square  matrix 


xy  - 

Z  . 

yy 


(5) 


which  may  be  regarded  as 
the  column  vector  ^  ^ 


the  mechanical  impedance  of  the  journal  bearing. 


e  is 


In  the  absence  of  other  constraints,  the  equation  of  motion  of  the  rotor  is 


m  s 


I  + 


0 


where  I  is  the  identity  matrix 


-1  0- 

-  0  1  - 


(6) 


Because  the  system  given  by  Eq.  (6)  is  homogeneous,  the  determinant  must  vanish, 


ms 


I  +  Z(8) 


(7) 


wheie 


m 


p  LD 
a 


The  root^  of  s  »  \  -f  iv  for  Eg.  (7)  represent  the  natural  exponential  constants 
of  the  system.  The  system  is  stable  or  unstable,  depending  on  whether  V  is  less 
or  greater  than  zero. 
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Expanding  Eq.  (7): 


(ras^  +  Z  )  (ms^  +  Z)-ZZ  »0 
XX  yy  xy  yx 


For  convenience,  introduce  the  definition 


z  =  u  +  iv  =  -  ras 
Z(s)  =  U  +  iV 


Then  Eq.  (8)  becomes 


z  (s),  s\  =  ^^2 


=  (Z  “  2 )(Z  -  z;  -  Z  Z  =  0 

XX  yy  xy  yx 


Separating  real  and  imaginary  parts  * 


W,  =  (U  -u)(U  -u)  -  (V  -v)(V  -v)-UU  +VV  »0 

1  XX  yy  XX  '  w  xv  vx  xv  vx 


XX  yy 


xy  yx  xy  yx 


W-  =  (U  -u)(V  -v)  +  (V  -v)(U  -u)-UV  -UV  »0 
2  XX  yy  '  xx  ''  yy  xy  yx  yx  xy 


At  the  state  of  neutral  stability 


s  *  iv 


u  s  m  V 
0  o 


V  =»  0 


By  analogy  with  a  single-degree-of- freedom  problem,  u  and  v  can  be  defined  as  the 
effective  spring  and  damping  of  the  system. 


*Alternately,  Equation  (11)  can  be  solved  for  z: 

z-i/z  +Z  +-\/(Z  -Z  )^  +  4Z  Z  > 

2  \  XX  yy  —  Y  '  XX  yy'  xy  yx  j 

t 

which,  however,  offers  no  advantage  over  Eq.  (12)  from  the  standpoint  of  ease  in 
computation. 
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"TT' 


Substituting  into  Eq.  (12): 

W,  (0,V  )  -  (U  -  m  V^)(U  -mV^)-VV  -UU  +VV  «0  (14) 

1  *  o'  '  XX  o  o''  yy  o  o'  xx  yy  xy  yx  xy  yx 

Wo(0,V  )  «  (U  -  m  V^)V  +  (U  -  m  V^)V  -  U  V  -  U  V  »  0  (15) 

2'  ’  o'  '  XX  o  o'  yy  '  yy  o  o  xx  xy  yx  yx  xy 


Solving  Eq.  (15)  for 


m 


UV  +UV  »UV  -UV 
XX  vv  vv  XX  XV  yx  yx  xy 

(V  +  V  ) 
o  '  XX  yy' 


(16) 


Substituting  into  Eq.  (14): 


[(U  -U  )V  +U  V  +U  V  lf-(U  -U  )V  +U  V  +U  V  I 
y  (0  V  )  »  L  XX  vv'  XX  XV  yx  vx  xyJL  '  xx  yy'  yy  xy  yx  yx  xy j 

^  °  (V  +  V  )^ 


XX  yy 


VV  -UU  +VV 
XX  yy  xy  yx  xy  yx 


(17) 


Since  all  elements  of  U  and  V  matrices  can  be  obtained  as  functions  of  v  by  solving 
the  pertinent  lubrication  equation,  Eq.  (17)  can  yield  roots  of  V^;  which  are  the 
characteristic  trequencies.  With  the  matrix  elements  evaluated  at  Eq.  (16) 
gives  the  critical  mass,  m^,  for  each  state  of  neutral  stability. 

For  a  slight  variation  from  the  state  of  neutral  stability,  writing 

8  ■  iv  +  68 
o 

m  ■  m  -f  8m 
o 


(18) 
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Sm  and  68  are  related  by  the  condition  that  6U  vanishes.  That  Is 


y  ds 


,  w  -  0 

Differentiating  Eq.  (9  ), 


(19) 


d  z  =  -  s  dm  -  2&sr  ds 


(20) 


Treating  z  as  an  analytic  function  given  by  Eq.  (11),  then  we  can  write 


dm  1  .dz. 

6»  -  0  '  '  8^ 


2m8) 


8W  »  0 


Because  Z  Is  directly  known  in  terms  of  V,  we  can  write 


dz 

ds 


8W  =  0 


i.  ^ 
I 


8W  a  0 


Then  at  s  =  iV, 


dm 

ds 


6W  a  0,  8  a  iv  iV 
o  o 


2m 


8W  -  0,  V 


V  1 
o  oJ 


'’^O 


8W  -  0,  V 


(21) 


Substituting  into  Eq.  (19),  for  slight  variation  from  the  state  of  neutral 
stability 


Iv  Eni 
o 


-  5 - j 

57  ■  V  > 


8W  -  0,  V 


Separating  real  and  imaginary  parts, 


5\  +  i8V 


{h  •  v^>)  +  <§7)^  J 


&m 


SW  -  0,  V 


(22) 


The  last  expression  is  an  explicit  stability  criterion  which  may  be  stated  as 
follows : 

With  an  inflnitessimal  mass  Increment  beyond  that  determined  from  the 
state  of  neutral  stability,  the  system  becomes  unstable  when  and  only 
when  evaluated  at  is  positive.  Conversely,  if  is 

negative  at  V^,  the  system  would  become  unstable  with  a  mass  decrement. 


To  find  ^u/^v  and  cVv/$V  ,  differentiate  Eq.  (11)  with  respect  to  V; 


h  +  [2'  -  “ 


Separating  real  and  imaginary  parts 


dv 


§7  • 


(23) 


2  2 

(q)  +  (r)" 


(24) 
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where 

q  -  &/(2q  - 

-  2u  -  (U  +  U  ) 

XX  yy 

r  +  Z^^)| 

=  2v  -  (V  +  V  ) 

'  XX  yy 

f  +Z  )  -  ItCZ  Z  -Z  Z  )\  (25) 

s  uU/^  dv'  XX  yy'  oy  xx  yy  xy  yx  /  ' 

»  i;^(U  +U  )  -  v|r(V  +V  )  -  l-CU  U  -V  V  -U  U  +  V  V  ) 
dv  XX  yy  ov  xx  yy  dV  xx  yy  xx  yy  xy  yx  xy  yx' 

-1  -  3«o|q|j<z^^4Zj^y)  -  57(z^^Zyy-z^yZy^)| 

*  -I?  <W’  ^  •  S7<'xxVVyy-“xyV''xyV> 


Slngle-Degree-of-Freedom  Systems 

While  above  derivations  were  performed  for  two-degree-of- freedom  systems,  the 
results  are  valid  for  slngle-degree-of- freedom  systems.  In  single^degree-of- 
freedom  systems;  e.g.,  parallel  thrust  plates,  and  a  non-spinning  externally 
pressurized  Journal  bearing, 


0 


(26) 


Equation  (6)  becomes  uncoupled.  Writing  Z  for  either  Z  or  Z  ,  Eq.  (7)  is 

XX  yy 

simply 


ms^  +  Z 


0 


(27) 
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The  state  of  neutral  stability  Is  then  given  by 

(28) 
(29) 


V(V^)  -  0 


A  slight  variation  from  the  state  of  neutral  stability  would  cause 


6s 


6\  +  i6v 


IV  ‘  8m 
o _ 

(3  2 


.  v*v 


r^V 


2  sv  2 


(30) 


U  and  V  clearly  have  similar  connotations  as  stiffness  and  damping  coefficients 
although  they  may  be  v  dependent  functions.  Above  results  have  previously  been 
given  in  connection  with  the  analysis  of  spiral -grooved  thrust  bearings  in 
Ref.  (19). 


Isotropic  Bearings 

For  those  bearings,  where  rotational  symmetry  prevails, 

(31) 

(32) 


XX 


yy 


.  Iff 


xy 


-  2i 


All  unloaded  cylindrical  Journal  bearings  are  in  this  category. 
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In  Ref. (24),  it  was  shown  that  reactions  of  an  Isotropic  bearing  caused  by  a  small 
periodic  motion  of  the  journal  are  linearly  related  to  those  due  to  a  small-orbit 
steady  whirl  of  the  journal.  Referring  to  Fig.  2,  the  bearing  reactions  of  the 
steady  whirl  problem  can  be  written  In  the  complex  notation  as  follows: 

(Fj^  +  IF^)  -  €  F  (A*)  (33) 

»  €  (f  +  ig) 

where  F(A*)  ■  f  -f  Ig  Is  a  complex,  analytical  function  of 

A*  -  A(l-2a')  (34) 

a*  "  +  V  respectively  for  the  forward  and  backward  whirl  cases.  The  impedance 
components  are  related  to  the  steady  %ihlrl  forces  In  the  following  manner  (see 
Ref.  (24)  for  details): 

Zjl  •  ■  V)  +  F(a'  -  -  V)]  (35) 

^ 

%diere  F^  Is  the  conjugate  of  F. 


The  characteristic  equation  Is  reduced  from  Bq.  (11)  to 


(*  -  Z,,)^  ■¥  Zj^  -  0 


II 


Expanding  and  making  use  of  Eqs.  (35)  and  (36): 


(37) 


-  ^F^(a'  -  V)  +  F(a’  -  -v)|  •  V)  F(a'  •  -v)  •  0 


(38) 
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At  the  state  of  neutral  stability 


V  •  V 


*  ■  u  « 


.v' 

O  O 


Then  Eq.  (38)  becomes: 


(39) 


From  the  imaginary  part  of  Eq.  (39), 


u 


2 

m  V 

o  o 


f(Vo)g(-V^)  -  f(-Vo)g(Vo) 

g(Vjj)  •  g(-v^) 


(40) 


Substituting  into  the  real  part  of  Eq>  (39): 

0  -  [f(v^)g(-v^)-f(-v^>g(v^)j^f(v^)g(-v^)-f(v^)g(v^) 

.2 

+  [»(v,)-t(-v„)J  [f(v^)  f(-v^)  +  g(v^)g(-v^)] 

•  [f(Vj^)l(-v,)-*(-v,)g(v^)]  [f(V,)g(V^)  -  «(-V,)g(-V^)] 

2 

+  [g(v,)  -  g(-v^)]  [j(v^)f(-v^)  +  g(v^)g(-v^)] 


The  solution  of  Eq.  (41)  is  determined  by  the  vanishing  of  either  of  three  factors. 
Vanishing  of  the  first  two  factors  require  that 


8  (±  Vq)  -  0  (42) 

Substituting  into  Eq.  (40) 

Vo  ■  <'‘3) 

Vanishing  of  the  third  factor  leads  to 

-  f(-v^)  -  g(v^)  -  g(-v^)  -  0  (44) 


This  condition  renders  Eq.  (40)  indeterminant.  Therefore,  m  has  to  be  found 

o 

by  other  means.  Substitute  Eq.  (44)  into  the  real  part  of  Eq.  (39),  it  is 
found 


There fore 


To  determine  whether  or  not  deviation  from  the  neutral  state  of  stability  in 
some  manner  would  cause  instability,  we  must  know  the  sign  of  dv/dv  at 
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Dlfferentlate  Eq.  (38)  with  respect  to  V: 


[2j  -  F^(v)  -  P(-v)]  ^  -  I  [f£(o'  -  V)  -  P(a'  .  -  V)] 

+  F(-v)  3^  -  F  (V)  ^.|  -  0 

a'  ■  V  I  a*  ■  -  V 


(45) 


Substituting  Eqs.  (42)  and  (43)  into  Eq.  (45), 


[f(t  '’„)-<(?  V^)  +  lg(+  v^)  ] 

-  f(±  V^)  fF^(a'  -  v)  -  F(oi'  -  -  V)] 


+  ''o>  55^ 


a'  -  V 


a' 


i8<+  55T  [f (a'  -!'')+  ig(a'  -tv)]  -  0 

[^(±  Vq)  -  f(+  v^)  +  ig(+  v^)]  [57  +  551  { f(a'  ■  ±  V)  +  ig(a' 


■  i  \>) ) 


-t 


V  •  V 


bgja'  ■  V) 
da^ 


V  •  V 


Separating  real  and  imaginary  parts 


(46 


V 

-  o 


■  •  i*. 


a*  •  -f  V 


(^7) 
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Uslng  these  results  In  Eq.  (22),  one  sees  that  the  stability  of  an  isotropic  bearing 
can  be  evaluated  in  terms  of  its  steady-whirl  properties.  Summarizing  Eqs.  (42), 
(43),  and  (47),  a  theorem  on  the  stability  of  an  isotropic  bearing  can  be  stated 
as  follows. 


Let  the  whirl  ratio  be  a'.  Let  the  radial  and  tangential  components  of  the  bearing 
reaction  due  to  the  whirl  motion  be  -  and  F^  respectively.  Let  €,  which  is 
assumed  to  be  small,  be  the  ratio  between  the  whirl  radius  and  the  radial  bearing 
clearance.  Then  the  state  of  neutral  stability  is  given  by 


0. 


t 

which  concerns  the  bearing  only  and  determines  0^  ;  and  also  by 


ffl 

o 


^-(a' ) 

R  o 

c(a;)2 


which  is  called  the  critical  mass. 


(48) 


(49) 


Let  the  rotor  mass  deviate  from  m^  by  a  sll^t  amount  da.  The  syatem  is  stable 
if  and  only  if 

00  >  0. 

Thus,  comparing  with  the  theorem  stated  below  Bq.  (22),  -F^  may  ba  regarded  as  the 
effective  damping  while  a’  now  assumes  the  role  of  v.  Note,  however,  a*  takes  on 
both  positive  and  negative  values  while  v  la  always  positive  by  definition. 


hf 
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PROCEDURE  OF  STABILITY  ANALYSIS 


The  essence  of  the  theorem,  and  its  simplified  versions,  derived  above,  will  be 
reiterated  in  the  form  of  a  procedure  of  stability  analysis.  The  dynamical 
model  considered  here  will  be  the  translational  mode  of  a  rigid  rotor  supported 
in  a  fixed  journal  bearing. 


The  gas  lubrication  equation  can  be  first  solved  for  the  equilibrium  position 
of  the  journal  designated  by  e^  as  shown  in  Fig.  1  along  with  other  parameters 
describing  the  operating  conditions,  e.g.,  A  and  L/D.  Then,  with  the  journal 
position  perturbed  by  the  periodic  motion 

«y)  (2) 

the  gas  lubrication  equation  can  again  be  solved  to  obtain  the  perturbed  forces, 
which  can  be  written  in  the  dimensionless  form: 


p  LD 


F  ) 

y 


(^x’  ®*P  (i'''*')! 


(50) 


The  numerical  values  of  (f  ,  f  )  depend  on  the  particular  geometry  as  well  as 

X  y 

the  equilibrium  operating  conditions  and  the  frequency  of  motion.  For  the 
present,  we  shall  merely  assume  that  the  dynamic  perturbation  forces  can  be 
found  for  the  entire  range  of  V  and  can  be  expressed  as 


(51) 
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The  elements  of  the  Impedance  matrix  are  complex  quantities 


Z  »  U  +  IV  etc, 

XX  XX  XX 


Case  1  €  ^  0.  Z  ^  G.  Z  ^0. 

-  o  xy  yx 


The  characteristic  frequencies  are  the  roots  of 


0  -  W^(v^) 


f(U  -U  )V  +  U  V  +  U  V  I  f-(U  -U  )V  +  U  V  +  u  V  1 
,,  r  XX  yy'  xx  xy  yx  yx  xyj  L  ^  xx  w  xv  vx  ^vx  xvl 


(V  +V  )' 
XX  yy 


VV  -UU  +VV 
XX  yy  xy  yx  xy  yx 


(17) 


The  critical  mass  is 


m 


UV  +UV  -UV  -HJ  V 
XX  yy  yy  xx  xy  vx  yx  xy 

V  ^  (V  +V  ) 
o  '  XX  yy 


(16) 


An  Inflnltesslmal  mass  increment  above  m^  would  constitute  an  unstable  system 
if 


2m  -  (U  +U  )]  fm  |^(V  +V  )  -  |;7(U  V  +V  U  -U  V  -V  U  )! 
o  o  XX  yy  J  I  o  o  oV'  XX  yy  dV'  xx  yy  xx  yy  xy  yx  xy  yx'J 


yy 


yy  XX  yy  xy  yx  xy  yx 


+  (V  +V  )  fm  1“  (U  +U  )  -  |^(U  U  -V  V  -U  U  +V  V  )] 
xx  yy'  I  o  o  oV  '  XX  yy  oV'  xx  yy  xx  yy  xy  yx  xy  yx'J 


>0 


at  V  »  V  ;  and  vice  versa, 
o 
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Case  3  Z  ■  Z  ■  Z//  ,  Z _ *  -  Z _ »  Zf 

-  XX  yy  "  *y  y*  ^ 

Generally,  this  corresponds  to  an  unloaded  Journal  bearing  possessing  rotational 
symmetry.  Zff  and  Zj_  are  directly  related  to  the  small  eccentricity,  steady- 
whirl,  radial  and  tangential  forces,  of  the  journal  bearing  with  whirl  ratio 
a'  «  +  V  (see  Fig.  2). 
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EXAMPLES 

In  Ref.  13,  the  stability  of  an  infinitely  long,  self-acting,  plain,  cylindrical 
Journal  bearing  vas  treated.  The  procedure  used  to  find  the  characteristic 
frequency  and  the  critical  mass  vas  essentially  the  same  as  Case  1  above.  The 
method  to  determine  which  side  of  the  neutral  state  is  unstable  was  different. 
Since  there  is  no  contradiction  in  fundamental  principle  between  the  two  methods, 
one  would  not  expect  the  final  results  to  be  different. 

Case  2  is  a  one-degree-of- freedom  problem.  The  corresponding  stability  criteria 
was  first  derived  in  Ref.  19,  which  also  gives  an  example. 

Two  examples  for  Case  3  will  be  considered  below. 


Example  1.  Unloaded,  Plain,  Cylindrical,  Self-Acting,  Journal  Bearing 


The  steady-whirl  results  of  this  bearing  was  given  in  Ref.  25.  They  are 
reproduced  in  FigS.  3  and  4.  The  tangentUl  hearing  force  vanishes  at  A*  =  0  or 

a'  =  0.5  for  all  A  and  L/D.  The  radial  bearing  force  is  also  always  zero 

o 

there.  Therefore,  the  critical  mass  is  zero.  Since,  at  A*  =  0, 


Since  all  finite  rotor  masses  are  larger  than  the  zero  critical  mass,  this 
bearing  is  always  unstable.  This  is  the  infamous  half- frequency  whirl. 
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Example  2.  Unloaded,  Herringbone -Grooved,  Cylindrical,  Self-Acting  .lournal  Bearing 


A  schematic  of  this  bearing  is  shown  in  Fig.  5  ,  its  steady-whirl  analysis  was 
given  in  Ref.  21.  Results  of  one  such  bearing  for  a  range  of  A  are  shown  in 
Figs.  6  and  7  .  The  value  of  a^,  at  which  the  tangential  force,  “F^,  vanishes, 
is  about  0.5  when  A  is  small;  it  Increases  with  A  and  is  about  1.7  at  A  »  50.  At 
a',  the  radial  force,  jF-,  has  some  positive  value,  which  also  Increases  with  A. 

O  K 

The  slope  of  — F^  at  is  always  negative.  Therefore,  masses  larger  than 


m  (A) 
o 


are  unstable,  while  smaller  masses  are  stable.  This  results  In  a  stability  map 
as  shown  in  Fig.  8. 
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CONCLUSIONS 

1.  The  response  of  a  gas  bearing  due  to  the  perturbation  of  a  periodic  motion 
can  be  used  to  determine  the  stability  map  of  the  bearing. 

2.  The  usual  concepts  of  stiffness  and  damping  can  be  generalized  to  include 
frequency  dependence  for  either  single-degree-of- freedom  or  coupled  systems. 

3.  The  condition  of  neutral  stability  exists  at  the  characteristic  frequency 
where  the  effective  system  damping  vanishes  and  when  the  mass  (<  Itlcal 
mass)  is  in  resonance  with  the  effective  system  stiffness. 

4.  Mass  increment  above  the  critical  mass  is  unstable  if  the  effective  system 
damping  Increases  with  frequency  and  conversely. 

5.  The  stability  map  of  unloaded  journal  bearings  can  be  constructed  from 
their  steady-whirl  results. 

6.  The  spectral  analysis  provides  a  satisfactory  means  to  separate  the  solution 
of  the  lubrication  equation  to  that  of  the  entire  dynamical  system. 
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NOMENCLATURE 


Bearing  Clearance 
Bearing  Diameter 

Displacement  of  Journal  Center  From  Bearing  Center 

Components  of  Rotor  Displacement  from  its  Equilibrium  Position 

Column  Vector  with  the  Components  f^  and  f^ 

Real  and  Imaginary  Parts  of  F 

Normalized  Complex  Moduli!  of  (F  ,  F  ) 

X  y 

Normalized  Complex  Bearing  Stiffness  for  Steady  Whirl 
(Normalization  Factor  is 

Conjugate  of  F 

Normalized  Radial  Bearing  Force,  ef 
Normalized  Tangential  Bearing  Force,  eg 
Components  of  the  Perturbed  Bearing  Force 
Idemfactor  or  Identity  Matrix 
Bearing  Length 

,  Nonaallzed  Rotor  Mass 

p^ld 

Rotor  Mass 
Ambient  Pressure 
Functions 
Bearing  Radius 

Notmalited  Complex  Exponential  Coefficient 
Time 

Real  Part  of  s 
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U 


U  ,  U  ,  etc. 
XX  xy 


V 


V  ,  etc. 

xy 


W 

W 

w 


1 

2 


z 


Z  ,  etc. 
xy 


(  ) 


0 


B(  ) 
a* 


e 


X 

A 

A* 


V 

T 

(0 


Real  part  of  Z 

Real  Parts  of  Z  ,  Z  .  etc. 

XX  xy 

Imaginary  Part  of  z 
Imaginary  Part  of  Z 
Imaginary  Parts  of  Z  ,  Z  •  etc. 

XX 

Characteristic  Determinant 

Real  Part  of  W 

Imaginary  Part  of  W 
2 

■ms  ,  Normalized  Effective  System  Impedance 

Normalized  Impedance  Hatrix  of  the  Bearing 

Components  of  the  Normalized  Impedance  Matrix 

Normalized  Parallel  or  Colinear  Impednace 

Normillzed  Perpendicular  or  Cross  Impednace 

Referring  to  the  Condition  of  Neutral  Stability  or  the 
Critical  Condition 

Inflnltesslmal  Deviation  of  (  )  from  Its  Critical  Value 

Whirl  Ratio,  Ratio  of  Angular  Whirl  Speed  to  u> 

Column  Vector  with  the  Components  c  and  c  > 

X  y 

e  /C,  e  /C 
X  y 

Real  Part  of  s 

Compressibility  Number 

P,  C 

A(l-2o') 

Viscosity 

lauginary  Part  of  a,  Normalised  Circular  Frequency  of  Oscillation 
Functions 

(ut.  Normalised  Time 
Rotor  Rotational  Speed 


)2+(c  ) 
X  '  y 
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